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1 Introduction 



Since time when quantum mechanics (QM) was born (in 1925-1926) two seem- 
ingly opposite approaches for description of atomic and subatomic physics were 
proposed respectively by Hciscnbcrg and Schrodingcr. Hciscnbcrg's approach is 
aimed at providing the affirmative answer to the fohowing question: Is combi- 
natorics of spectra (of observables) provides sufficient information about mi- 
croscopic system so that dynamics of such a system can be described in terms of 
known macroscopic concepts? Schrodinger's approach is exactly opposite and is 
aimed at providing the aSirmative answer to the following question: Using some 
plausible mathematical arguments is it possible to find an equation which under 
some prescribed restrictions will reproduce the spectra of observables? Although 
it is widely believed that both approaches are equivalent, already Dirac in his 
lectures on quantum field theory [1] noticed (without much elaboration) that 
Schrodinger's description of QM contains a lot of "dead wood" which can be 
safely disposed altogether. According to Dirac "Heisenberg's picture of QM is 
good because Hciscnbcrg's equations of motion make sense" . To our knowledge, 
Dirac's comments were completely ignored, perhaps, because he had not pro- 
vided enough evidence making Heisenberg's description of QM superior to that 
of Schrodinger's. In recent papers [2,3] we found examples supporting Dirac's 
claims. From the point of view of combinatorics, there is not much difference 
in description of QM, quantum field theory and string theory. Therefore, in 
this paper wc choose Hciscnbcrg's point of view on string theory using results 
of our recent works in which we re analyzed the existing treatments connecting 
Veneziano (and Veneziano-like) amplitudes with the respective string-theoretic 
models. As result, wc were able to find new tachyon-free models reproducing 
Veneziano (and Veneziano-like) amplitudes. In this work results of our papers 
[4-7] which will be respectively called as Part I, Part II, Part III and Part IV 
are further developed to bring them in correspondence with those developed 
by other authors. Without any changes in the already existing formalism, we 
were able to connect our results with an impressive number of string-theoretic 
models. 

As is well known, all information in high energy physics is obtainable through 
proper interpretation of the scattering data. It is believed that for suSiciently 
high energies such data arc well described by the phcnomcnological Regge theory 
and can be conveniently summarized with help of Chew-Frauthchi plots relating 
masses to spins (angular momenta), e.g. see book by Collins[8]. Veneziano 
amplitudes are by design Regge-behaving. Both Regge theory and Veneziano 
amplitudes emerged before major developments in QCD took place in 70ies. 
Once these developments took place, naturally, it was of interest to recover the 
Regge theory from QCD. Even though there are many ways for doing so, to 
our knowledge, the problem is still not solved completely. This is so for the 
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following reasons. 

Although the amount of data obtained by perturbative treatments of QCD is 
quite impressive, e.g. read Ref.[9], these results are not helpful for establishing 
the Regge-type behavior of QCD. Since such a behavior can be easily established 
with help of variety of string models, the task lies is connecting these models 
with QCD. Evidently, such a task is also equally applicable to models reproduced 
in this work. 

The connections between either QCD and spin chains or between strings 
and spin chains were already discussed in literature for quite some time. Re- 
cent paper by Dorey [10] contains may references listing these earlier results. 
Subsequently, they had been replaced by those whose methods are based on 
ADS/CFT correspondence. Prom the point of view of this correspondence, the 
connections between strings and QCD also can be made through spin chains 
as it is demonstrated in the seminal paper by Gubser, Klebanov and Polyakov, 
Ref.[ll]. Their ideas were developed in great detail in the paper by Minahan 
and Zarembo, Ref.[12]. The spectrum of anomalous dimensions of operators 
in the A/" = 4, N—^ oo supersymmetric Yang-Mills (Y-M) model (described in 
terms of the excitation spectrum of the spin chain model) is related to the 
string spectrum describing hadron masses. To connect these facts with devel- 
opments in this paper we mention papers by Kruczenski [13], and Cotrone et 
a/.[14]. In both papers spin 1/2 XXX Heisenberg chain was used for description 
of the excitation spectrum. Furthermore, in the paper by Cotrone et al the 
explicit connection with the hadron mass spectrum was made. Both papers 
use ADS/CFT correspondence in their derivations. In this work, we reobtain 
these spin chain results using combinatorial arguments following Heisenberg's 
philosophy discussed in some detail in our recent work [2]. In view of this, our 
derivation is not relying on use of ADS/CFT methods. 

It should be noted though that the ADS/CFT is not the only method of 
connecting QCD with strings. In 1981 't Hooft suggested [15] to reduce the non 
Abelian QCD to Abelian Ginzburg-Landau (G-L) type theory. The rationale for 
such an Abelian reduction can be traced back to the work by Nambu [16]. In his 
work Nambu superimposed G-L theory with the theory of Dirac monopoles to 
demonstrate quark confinement for mesons. Incidentally, Veneziano amplitudes 
are suited the most for describing meson resonances, e.g. see Ref.[8]. These 
are made of just two quarks: quark and antiquark. Thus, if the existence of 
Abelian reduction would be considered as proven, this then would be equivalent 
to the proof of quark confinement. Recent numerical studies have provided con- 
vincing evidence supporting the idea of quark confinement through monopole 
condensation, e.g. see Refs. [17,18]. Since publication of 't Hooft's paper many 
theoretical advancements were made, most notably by Cho, Ref.s[19,20], and 
Kondo, Ref.s [21-23], whose work was motivated by that by Faddeev and his 
group. Results of this group are summarized in the recent review by Faddeev, 
Ref . [24] . From this reference it follows that most of efforts to date were spent on 
description of the massless version of QCD. The excitation spectrum of solitonic 
knotted-like structures (admitting interpretation in terms of closed strings) pro- 
vides the spectrum of glueball masses as demonstrated in [22]. In the recent 
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paper [25] Auckly, Kapitanski and Speight demonstrated how Skyrme model 
can be obtained from Faddeev model. Since Skyrme model was used for a long 
time for description of the baryon spectra, e.g. read Ref.[26], and since already 
Nambu recognized usefulness of the Abclian reduction for description of the 
meson spectra, it follows that the Abelian reduction method is also capable of 
providing sufficient information about QCD in the strong coupling regime. 

Since our derivation is not relying on use of ADS/CFT methods, naturally, 
the question arises: Can Heisenberg-style analysis of this work be useful in de- 
ciding which of the methods: ADS/CFT or Abelian reduction is more likely 
realized in Nature? In our companion paper [27] we provide evidence of superi- 
ority of the Abelian reduction over ADS/CFT methods. To avoid duplications, 
in this work we only develop connections between Veneziano amplitudes, spin 
chains and various string models without making attempts to relate explicitly 
these models with QCD. 

The rest of this paper is organized as follows. In Sections 2-3 we provide 
needed physical motivations and background. For this purpose, in Section 
2 we replaced mathematically sophisticated derivation of Veneziano partition 
function developed in [4-6] by a considerably simpler combinatorial derivation 
of this fimction analogous to that discussed in [7]. As a by-product of this effort, 
we were able to uncover connections with spin chains already at this stage of 
our investigation. To strengthen this connection, in Section 3 we demonstrate 
that the obtained Veneziano partition function coincides with that for the 
Polychronakos-Frahm (P-F) partition function for the ferromagnetic spin chain 
model [28]. In the limit of infinitely long chains such a model coincides with 
spin 1/2 XXX Heiscnberg spin chain. This is explained in Section 4. In the same 
section we discuss different paths aimed at establishing links between the P-F 
spin chain and variety of string-theoretic models, including the most recent ones. 
This is achieved by mapping the combinatorial and analytical properties of the 
P-F spin chains into analogous properties of spin chains used for description 
of the stochastic process known in literature as asymptotic simple exclusion 
process (ASEP). To make our presentation self-contained. Appendix contains 
some basic information on ASEP sufficient for understanding the results of the 
main text. In addition, in the main text, we provide some information on the 
Kardar-Parisi-Zhang (KPZ) and Edwards- Wilkinson (EW) equations which are 
just different well defined macroscopic limits of the microscopic ASEP equations. 
This is done with purpose of reproducing a variety of already known string- 
theoretic models, including the most recent ones, and to relate them to each 
other. Section 5 contains some discussion connecting the obtained results with 
that presented in our companion paper [27]. 
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2 Combinatorics of Veneziano amplitudes and 
spin chains: qualitative considerations 

In Part I, we noticed that the Veneziano condition for the 4-particle amphtude 
given by 

a{s) + a{t) + a{u) ^ -1, (1) 

where a(s), a{t), a{u) € Z, can be rewritten in more mathematicaUy suggestive 
form. To this purpose, following [29], we need to consider additional homogenous 
equation of the type 

a{s)m + a{t)n + a{u)l + k ■ 1 ^ Q (2) 
with m, n, I, k being some integers. By adding this equation to (1) we obtain 

a{s)'m + a{t)n + a{u)l = k (3a) 

or, equivalently, 

"■1 + "-2 + "-3 = N , (3b) 

where all entries hy design are nonnegative integers. For the multiparticle case 
this equation should be replaced by 

no + --- + nk = N (4) 

so that combinatorially the task lies in Ending all nonnegative integer combi- 
nations of nQ,...,nk producing (4). It should be noted that such a task makes 
sense as long as N is assigned. But the actual value of N is not fixed and, hence, 
can be chosen quite arbitrarily. Eq.(l) is a simple statement about the energy 
-momentum conservation. Although the numerical entries in this equation can 
be changed as we just explained, the actual physical values can be subsequently 
re obtained by the appropriate coordinate shift. The arbitrariness in selecting 
N reflects a kind of gauge freedom. As in other gauge theories, we may try to 
fix the gauge by using some physical considerations. These include, for example, 
an observation made in Part I that the 4 particle amplitude is zero if any two 
entries into (1) are the same. This fact prompts us to arrange the entries in 
(3b) in accordance with their magnitude, i.e. ni > n2 > n^. More generally, we 
can write: > rii > • • • > n^, > 10. 

Provided that (4) holds, we shall call such a sequence a partition and shall 
denote it as n =(no, ...^Uf^). If n is a partition of iV, then we shall write n\- N. 
It is well known [30,31] that there is one- to -one correspondence between the 
Young diagrams and partitions. We would like to use this fact in order to 
design a partition function capable of reproducing Veneziano (and Veneziano- 
like) amplitudes. Clearly, such a partition function should also make physical 
sense. Hence, we would like to provide some qualitative arguments aimed at 

^The last inequality; nj. > 1, is chosen only for the sake of comparison with the existing 
literature conventions, e.g. see Ref. [30]. 
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convincing our readers that such a partition function does exist and is physically 
sensible. 

We begin with observation that there is one- to- one correspondence between 
the Young tableaux and directed random walk^. It is useful to recall details of 
this correspondence now. To this purpose we need to consider a square lattice 
and to place on it the Young diagram associated with some particular partition. 

Let us choose some h x m rectangle^ so that the Young diagram occupies 
the left part of this rectangle. We choose the upper left vertex of the rectangle 
as the origin of the xy coordinate system whose y axis (South direction) is 
directed downwards and x axis is directed Eastwards. Then, the South-East 
boundary of the Young diagram can be interpreted as directed (that is without 
self intersections) random walk which begins at (0, — m) and ends at {n,0). 
Evidently, such a walk completely determines the diagram. The walk can be 
described by a sequence of O's and I's. Say, for the x— step move and 1 for the 
y— step move. The totality Af of Young diagrams which can be placed into such 
a rectangle is in one-to-one correspondence with the number of arrangements 
of O's and I's whose total number is m -|- n. Recalling the Fermi statistics, the 
number Af can be easily calculated and is given hy Af ~ (m + n)!/m!n|l. It can 
be represented in two equivalent ways 



(m -|- n)!/m!n! 



(n + l)(n + 2) • • • (n + m) _ / n + m 

m\ ~ \ m 

{m + l){in + 2) ■ ■ ■ {n + m) _ f m + n 



(5) 



Let now p{N; fc, m) be the number of partitions of N into < k nonnegative 
parts, each not larger than ni. Consider the generating function of the following 
type 



J-{k, m I q) 



s 



N=0 



p{N; fc, m)q 



N 



(6) 



where the upper limit S will be determined shortly below. It is shown in 

k + m _ k + m 
m ~ k 



Refs. [30,31] that T{k, m \ q) 



where, for instance, 



k -\- m 
m 



0. From this result it should be clear that the expression 



k + m 
m 



k + m 
m 



•^Furthermore, it is possible to map bijectively such type of random walk back into Young 
diagram with only two rows, e.g. read [32], page 5. This allows us to make a connection with 
spin chains at once. In this work we are not going to use this route to spin chains in view of 
the simplicity of alternative approaches discussed in this section. 

^Parameters h and m will be specified shortly below. 

^We have suppressed the tildas for n and m in this expression since these parameters are 
going to be redefined below anyway. 

^On page 15 of the book by Stanley [31], one can find that the number of solutions -/V(n, k) 

n + k-1 
k-1 



in positive integers to j/i + ... + i/^ = n + fc is given by 



while the number of 



9=1 



6 



is the g— analog of the binomial coefficient 

_ ((?«-l)((7«-i-l) 



In literature [30,31] this 



k + m 
m 

analog is known as the Gaussian coefficient. Explicitly, it is defined as 



h-1 



l)...(q-l) 



(7) 



for some nonegative integers a and b. From this definition we anticipate that 
the sum defining generating function !F{k^ m \ q) in (6) should have only finite 
number of terms. Eq.(7) allows easy determination of the upper limit S in the 
sum (6). It is given by km. This is just the area of the k x m rectangle. In 
view of the definition of p{N; k,m), the number m — N ~ k. Using this fact 

N 



(6) can be rewritten as: T{N^ k \ q) 



k 



.This expression happens to be 



the Poincare' polynomial for the Grassmannian Gr(m, k) of the complex vector 
space C^of dimension N as can be seen from page 292 of the book by Bott 
and Tu, [33j3. From this (topological) point of view the numerical coefficients, 
i.e. p{N; fc, to), in the q expansion of (6) should be interpreted as Betti numbers 
of this Grassmannian. They can be determined recursively using the following 
property of the Gaussian coefficients [31], page 26, 



(8) 



= 1. We refer our readers to Part II for 



' n + 1 ' 




n 




n 


k + 1 




k + 1 




k 




q 




q 





and taking into account that 



n 




mathematical proof that T{N, k\ q) is indeed the Poincare' polynomial for the 
complex Grassmannian. With this fact proven, we notice that, due to relation 
TO = — fc, it is sometimes more convenient for us to use the parameters to 
and k rather than N and k. With such a replacement we obtain: 



!F{k^ TO 



k 

i=l 



k + m 
k 



,m+l 



1) 



l)...(g-l) 



(9) 



This result is of central importance for this work. In Part II a considerably more 
sophisticated mathematical apparatus was used to obtain it (e.g. see equation 
(6.10) of this reference and arguments leading to it). 

In the Hmit q ^ 1 Eq.(9) reduces to A/" as required. To make connections 
with results known in physics literature we need to re scale q's in (9), e.g. let 



solutions in nonnegative integers to xi 



Careful reading of Page 



n + k 
k 

15 indicates however that the last number refers to solution in nonnegative integers of the 
equation xq + ... + Xj, = n. This fact was used essentially in (1.21) of Part I. 

^To make a comparison it is sufficient to replace parameters and n in Bott and Tu book 
by q and A^. 
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q = t'i . Substitution of such an expression back into (9) and taking the hmit 
t — *■ 1 again produces J\f in view of (5). This time, however, we can accomphsh 
more. By noticing that in (4) the actual value of deliberately is not yet fixed 
and taking into account that m = N — k we can fix N by fixing m. Specifically, 
we would like to choose m — 1 ■ 2 ■ 3 ■ ■ ■ k and with such a choice we would like 
to consider a particular term in the product (9), e.g. 

f - t^+'r 

Si^) = . (10) 

In view of our "gauge fixing" the ratio m/i is a positive integer by design. This 
means that we are having a geometric progression. Indeed, if we rescale t again 
: t ~* t^, we then obtain: 

5(^) = l + ^2 + ... + ^2A (jj) 

with rh — Written in such a form the above sum is just the Poincare' poly- 
nomial for the complex projective space CP™. This can be seen by comparing 
pages 177 and 269 of the book by Bott and Tu [33]. Hence, at least for some 
m's, the Poincare' polynomial for the Grassmannian in just the product of the 
Poincare' polynomials for the complex projective spaces of known dimensional- 
ities. For m just chosen, in the limit i — > 1, we reobtain back the number J\f 
as required. This physically motivating process of gauge fixing just described 
can be replaced by more rigorous mathematical arguments. The recursion re- 
lation (8) indicates that this is possible. The mathematical details leading to 
factorization which we just described can be found, for instance, in the Ch-3 of 
lecture notes by Schwartz [34]. The relevant physics emerges by noticing that 
the partition function Z{J) for the particle with spin J is given by [35] 

Z{J) = trie-f^"^"^) = ef-' + e'''^-^-^^ + ■ ■ ■ + e-"-' 

= 6^=^(1 + e-"-Fe-2<= + ... + e-2^''), (12) 

where c is known constant. Evidently, up to a constant, Z[J) ~ S{i). Since 
mathematically the result (12) is the Weyl character formula, this observation 
brings the classical group theory into our discussion. More importantly, because 
the partition function for the particle with spin J can be written in the language 
of N=2 supersymmetric quantum mechanical modeH, as demonstrated by Stone 
[35] and others [36] , the connection between the supersymmetry and the classical 
group theory is evident. It was developed in Parts III and IV. 

In view of arguments presented above, the Poincare' polynomial for the 
Grassmannian can be interpreted as a partition function for some kind of a spin 
chain made of apparently independent spins of various magnitudetQ. These 
qualitative arguments we would like to make more mathematically and physi- 
cally rigorous. The first step towards this goal is made in the next section. 

*We hope that no confusion is made about the meaning of N in the present case. 

^In such a context it can be vaguely considered as a variation on the theme of the Polyakov 
rigid string (Grassmann a model, Ref. [37], pages 283-287), except that now it is exactly 
solvable in the qualitative context just described and, below, in the mathematically rigorous 
context. 



8 



3 Connection with the Polychronakos-Frahm spin 
chain model 



The Polychronakos-Frahm (P-F) spin chain model was originally proposed by 
Polychronakos and described in detail in [28]. Frahm [38] motivated by the 
results of Polychronakos made additional progress in elucidating the spectrum 
and thermodynamic properties of this model so that it had become known as 
the P-F model. Subsequently, many other researchers have contributed to our 
understanding of this exactly integrable spin chain model. Since this paper is 
not a review, we shall quote only works on P-F model which are of immediate 
relevance. 

Following [28] , we begin with some description of the P-F model. Let af (a = 
1, 2, — 1) be SU{n) spin operator of i-th particle and let the operator cr.y 
be responsible for a spin exchange between particles i and j, i.e. 

a,, = l+Y.<<^';- (13) 

In terms of these definitions, the Calogero-type model Hamiltonian can be writ- 
ten as [39,40] 

w = + + (14) 

^ i i<j (Xi - Xj) 

where I is some parameter. The P-F model is obtained from the above model 
in the limit I —> ±oo . Upon proper rescaling of Ti. in (14), in this limit one 
obtains 

nv-:F - -signil) E (15) 

where the coordinate Xi minimizes the potential for the rescaled Calogero models, 
that is 



^'^^ = E 73- (16) 

i<j (Xi - Xj) 

It should be noted that TL-p^jr is well defined without such a minimization, 
that is for arbitrary real parameters Xi. This fact is explained in detail in our 
companion paper [27]. In the large / limit the spectrum of TL is decomposable as 

E = Ec+ lEv-r. (17) 

where Ec is the spectrum of the spinless Calogero model while E-p^jr is the 
spectrum of the P-F model. In view of such a decomposition, the partition 
function for the Hamiltonian Ti. at temperature T can be written as a product: 
Z-h{T) =Zc{T)Z-p^jr(T / 1) . From here, one formally obtains the result: 

Zp_^(T) = hm ^iM. (18) 



^''The Calogero model is obtainable from the Hamiltonian (14) if one replaces the spin 
exchange operator aij by 1. Since we are interested in the large I limit, one can replace the 
factor l{l — I) by P in the interaction term. 
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It implies that the spectrum of the P-F spin chain can be obtained if both the 
total and the Calogero partition functions can be calculated. In [28] Polychron- 
akos argued that Zc(T) is essentially a partition function of noninteracting 
harmonic oscillators. Thus, we obtain 

^ 1 1 

Zc{N;T) = n 1 9 = exp(-/3),/3 = {ksT)-^ ■ (19) 

i=i ^ — q 

Furthermore, the partition function Zfi{T) according to Polychronakos can be 
obtained using Zc{N;T) as follows. Consider the grand partition function of 
the type 

^ N — I ^ ry IT .rp\„,L 



S - ^ ^„(iV;T)2/^ = E Zc{L-T)y^ (20) 
where n is the number of flavorJ^. Using this definition we obtain 

n 

Z„(iV;T)= V^^c{h;T). (21) 

Next, Polychronakos identifies Z„(iV;T) with Zh(T). Then, with help of (18) 
the partition function Z-p^jr(T) is obtained straightforwardly as 

N 

n(i-9') 

Zr-AN;T)= Y ■ (22) 

i—1 r—1 

Consider this result for a special case: n — 2. It is convenient to evaluate the 
ratio first before calculating the sum. Thus, we obtain 



2 fe^ {l-q)---{l~q'^){l~q)---(l-q^-'') 

n n(i-'7'-) 

i—l r—1 

(23) 

where the Poincare' polynomial J-'{N, k \ q) for the Grassmanian of the complex 
vector space of dimension N was obtained in the previous section. Indeed 
(23) can be trivially brought into the same form as that given in our Eq.(9) 
using the relation m + k = N. To bring Eq.(9) in correspondence with equation 
(4.1) of Polychronakos [28], we use the second equality given in Eq.(9) in which 
we make a substitution: m = N — k. After this replacement, Eq.(22) acquires 
the form 

N k ^ N-i+l 

Z^_AN,T) = Y\{ 7_ . (24) 

k=0 'i=0 ^ 



-^-"^That is n the same number as n in SU(n). 
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coinciding with Eq.(4.1) by Polychronakos. This equation corresponds to the 
ferromagnetic version of the P-F spin chain model. To obtain the antiferromag- 
netic version of the model requires us only to replace q by in (24) and to 
multiply the whole r.h.s. by some known power of q. Since this factor will not 
affect thermodynamics, following Frahm [38], we shall ignore it. As result, we 
obtain 

^(iV; T) = ^ qW2-^r n 7 . , (25) 

fe=0 i=0 

in accord with Frahm's Eq.(21). This result is analyzed further in the next 
section playing the central role for this paper. 

4 Connections with WZNW model and XXX 
s=l/2 Heisenberg antiferromagnetic spin chain 

4.1 General remarks 



To establish these connections we follow work by Hikami [41] . For this purpose, 
we introduce the notation 



(9)n=n(i-'?') 



(26) 



allowing us to rewrite (25) in the equivalent form 



k=0 



i=0 



q' 



fc=0 



(N/2-kf 



N 



N-k 



(27) 

Consider now the limiting case (A'' — > oo) of the obtained expression. For this 
purpose we need to take into account that 



lim 

N^oo 



N 
k 



1 

Wk' 



(28) 



To use this asymptotic result in (27) it is convenient to consider separately 
the cases of N being even and odd. For instance, if N is even, we can write: 
N = 2m. In such a case we can introduce new summation variables: I — k ~ m 
and/or I = m — k. Then, in the limit N ^ oo (that is m ^ oo) we obtain 
asymptotically 



Z;i^(oo;r) = ^ ^ q^ . (29a) 
in accord with [41]. Analogously, if N = 2m + 1, we obtain instead 

-. oo 
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According to Melzer [42] and Kedem, McCoy and Melzer [43], the obtained 
partition functions coincide with the Virasoro characters for SUi(2) WZNW 
model describing the conformal hmit of the XXX (s=l/2) antiferromagnetic 
spin chain [44]. Even though equations (29a) and (29b) provide the final result, 
they do not reveal their physical content. This task was accomplished in part 
in the same papers where connection with the excitation spectrum of the XXX 
antiferromagnetic chain was made. To avoid repetitions, below we arrive 
at these conclusions using different arguments. By doing so many new and 
unexpected connections with other stochastic models will be uncovered. 

4.2 Method of generating functions and q-deformed har- 
monic oscillator 

We begin with definitions. In view of equations (9), (24) and (27), we would 
like to introduce the Galois number Gn via 

(30) 

This number can be calculated recursively as it was shown by Goldman and 
Rota [45] with the result 

Gn+1 = 2Gn + {q^ - 1) Gn-1. (31) 

Alternative proof was given by Kac and Cheung [46]. To calculate Gat we 
have to take into account that Gq ~ 1 and Gi = 2. These results can be used 
as a reference when one attempts to calculate the related Rogers-Szego (R-S) 
polynomial Hi^{t) defined as [47] 

t'' (32) 

so that Hpf{l) — GjvEl. Using Ref.[46] once again, we find that 11^(1) obeys 
the following recursion relation 

H^+iit) = {l + t)HN{t) + {q^ - 1) tH^-iit) (33) 

which for t = 1 coincides with (31) as required. The above recursion relation is 
supplemented with initial conditions. These are : Hq — 1 and Hi = 1 + t. 

At this point we would like to remind our readers that for t = 1 according to 
(24), (30) and (32) we obtain: Z^_jr{N;T) = Gn- Hence, by calculating i?Ar(i) 
we also shall obtain the partition function for the P-F chain. 

To proceed with such calculations, we follow Ref.[48]. In particular, we 
consider first the auxiliary recursion relation for Hermite polynomials: 

Hn+iix) = 2xHn{x) - 2nHn-i{x) (34a) 

-'^For brevity, unless needed explicitly, we shall suppress the argument q in Hi^{t;q). 



Gat = ^ 



J, n 



N 
k 



HM{t-q) :=^ 
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supplemented by the differential relation 



—Hn{x) = 2nHn-i{x) 
which, in view of (34a), can be conveniently rewritten as 



(34b) 



Hn+l{x) = {2x-—)Hn{x). 



(34c) 



This observation prompts us to introduce the raising operator R = 2x ^ so 

ax 

that we obtain: 



R"Ho{x)=H,,{x). 



(35) 



The lowering operator can be now easily obtained again using (34). Indeed, we 
obtain 

~Hn{x) = LHn{x) = nHn-i{x) (36) 

so that [L, i?] = 1 as required. Based on this, the number operator N can be 
obtained as = RL so that NHn{x) = nHn{x) or, explicitly, using provided 
definitions, we obtain: 



i^-2xl + 2n)H.ix)=0. 



(37) 

1 > and , 



Evidently, we can write: R \ n >—\ n + 1 >, L \ n >=\ n 
< m I n n\S„in as usual. 

Wc would like now to transfer all these results to our main object of interest- 
the recursion relation (33). To this purpose, wc introduce the difference operator 
A via 

AHN{t) := H^it) - HN{qt). (38) 
Using (32) we obtain now 

AHN{t) = (l-g^)tifjv-i(i), 
where we took into account that 



(39) 



" N ' 




N 


k 


Q 


_ N-k _ 



Using this result in (33) we obtain at once 

HN+i{t) = [{l + t)-A]HNit). 



(40) 



(41) 



This, again, can be looked upon as a definition of a raising operator so that we 
can formally rewrite (41) as 



nHN{t) = HN+i{t). 



(42) 
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The lowering operator can be defined now as 

£ := -A (43) 

X 

so that 

CHN{t) = {l-q'')HN-iit). (44) 
The action of the number operator Af = TZC is now straightforward, i.e. 

J^HN{t)^{l-q'')HN{t). (45) 

Following Kac and Cheung [32] we introduce the (j— derivative via 

(46) 

By combining this result with (38) we obtain, 

This allows us to rewrite the raising and lowering operators in terms of g— derivatives. 
Specifically, we obtain 

n :=(1 + i) - (1 - q)tDg, (48) 

and 

C := Dg. (49) 

While for the raising operator rewritten in such a way Eq.(42) still holds, for 
the lowering operator £ we now obtain: 

CHN{t) = -Y^HN-i{t) = [N]HM-i{t). (50) 

The number operator J\fq is acting in this case as 

NqHN{t) = [N]HN{t). (51) 

Wc would like to connect these results with those available in literature on 
g— deformed harmonic oscillator. Following Chaichan et al [49], we notice that 
the undeformed oscillator algebra is given in terms of the following commutation 
relations 

aa^ — a^a = 1 (52a) 
[iV, a] = -a (52b) 

and 

[N,a+]=a+. (52c) 

In these relations it is not assumed a priori that A'' = a+a and, therefore, 
this algebra is formally different from the traditionally used [a, a+] = 1 for 
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the harmonic oscillator. This observation allows us to introduce the central 
element Z = N — a^a which is zero for the standard oscillator algebra. The 
deformed oscillator algebra can be obtained now using Eq.s(52) in which one 
should replace (52a) by 

aa^ — qa^ a = 1. (52d) 

as it is done in [50]. Consider now the combination K := ClZ-qN acting on Hm 
using previously introduced definitions. Simple calculation produces an operator 
identity: CTZ — qAf =1, so that we can formally make a provisional identification 
: £ and TZ ^a"*". To proceed, we need to demonstrate that with such an 
identification equations (52 b,c) hold as well. For this to happen, we should 
properly normalize our wave function in accord with known procedure for the 

harmonic oscillator where we have to use | n >= —= (a'*')" | >. In the present 

vn! 

case, we have to use I N >= — r= (TZ)" I > as the basis wavefunction while 

making an identification: | >= Ho{t). The eigenvalue equation (51), when 
written explicitly, acquires the following form: 

[tDl - Y^Dg + ^^]HN{t) = 0. (53) 

4.3 The limit g — > 1^ and emergence of the Stieltjes-Wigert 
polynomials 



Obtained results need further refinements for the following reasons. Although 
the recursion relations (33), (34) look similar, in the limit q ^ 1^ equation (33) 
is not transformed into (34). Accordingly, (53) is not converted into equation 
for Hermite polynomials known for harmonic oscillator. Fortunately, the situ- 
ation can be repaired in view of recent paper by Karabulut [51] who spotted 
and corrected some error in the influential earlier paper by Macfarlane [52]. 
Following the logic of these papers, we define the translation operator T(s) as 
T(s) := e^v^ . Using this definition, the creation and annihilation a operators 
are defined as follows 

at = J. [(7^+3 - {s)]T^ (s) (54a) 

V 1 ^ 9 

where T^{s) = and, accordingly, 

a = -^Ti{s)[q^+i - Ti{s)]. (54b) 
Under such conditions, the inner product is defined in the standard way, that is 

oo 

if, 9)= I rix)g{x)dx (55) 
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so that {q^)^ = and (d/dx)^ = —{d/dx) thus making the operator to 
be a conjugate of a in a usual way. The creation-annihilation operators just 
defined satisfy commutation relation (52d). At the same time, the combination 
a^a while acting on the wave functions (to be defined below) produces 
equation similar to (51), that is 



Furthermore, it can be shown, that 

a*n = \/K'^n-i and a^*„ = ^/Xn+l'^'n+l 



(56) 



(57) 



in accord with previously obtained results. Next, we would like to obtain the 
wave function \l/„ explicitly. To this purpose we start with the ground state 
a^'o = and use (54b) to get (for s — 1/2 the following result 



*0 {x+^) = qi+-^o (x). 



(58) 



Let w{x) be some yet unknown function. Then, it is appropriate to look for 
solution of (58) in the form 



^'o (x) — const ■ w{x)q^ 



(59a) 



provided that the function w{x) is periodic: w{x) — 'w{x+l/2). The normalized 
ground state function acquires then the following look 



where the constant a-w is given by 



(59b) 



/ oo 








w{x) 




\— oo 







(59c) 



Using this result, can be constructed in a standard way through use of the 
raising operators. There is, however, a faster way to obtain the desired result. 
To this purpose, in view of (59b), suppose that (x) can be decomposed as 
follows 



(60) 



Viq^q)n k=0 

where {q,q)n = (1 — <7)(1 — q^) • • • (1 — q") and CJ}{q) is to be determined. 
To do so, by applying the operator / yj A„+i to (60) and taking into account 
that T^2 w(x) — w(x — 1/2) = w{x) (in view of the periodicity of ii;(x)) we 
end up with the recursion relation for CJ!{q): 



Cl+\q)=q'^Cl{q) + Cl_M- 



(61) 



^The rationale for choosing s = 1/2 is explained in the same reference. 
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This relation should be compared with that given by (8). Andrews [47], page 
35, demonstrated that (8) and (61) are equivalent. That is, 



CUq)= I ■ (62) 

J q 

This result implies that, indeed, up to a constant, the obtained wavefunction 
should be related to the Rogers-Szego polynomial. The obtained relation is 
nontrivial nevertheless. We would like to discuss this nontriviality in some 
detail now. 

Following [51,53], let q = e"'^ , where c is some nonegative number. Intro- 
duce the distributed Gaussian polynomials via 

oo 

<i>„(x) = ^C•,"(g)(-l)'=g-'=/2g(-'=)^ (63) 

fe=0 

These polynomials satisfy the following orthogonality relation: 



<^n{x)<^m{x)dx = \\<^n{x)f (64) 



witlilfl 

\\'^n{x)\\ = [-^f q-"^ Vi^n. (65) 

This result calls for a change in normalization of $„(x), i.e.,(/)„(x) = ||^"|^||| . 
Under such conditions 4>n{x) coincides with (x), provided that w{x) — 1. 
Introduce new variable u = q~^^, and consider a shift: $„(x) ^n{x — s). 
Using (63), we can write 

$„(x-s) ^ exp{~ {\nuf / {-4:\nq)}P„{u; s), (66a) 

where 

oo 

P„(u; s)^Yl Cfe"(<z)(-l)'=<?-'=/'<z(^+'=)'M^ (66b) 

fe=0 

The orthogonality relation (64) is converted then into 

oo 

j dww2«~iexp{-(lnu)V(-41ng)}P„(u;s)P„(u;s) = ,5™„. (67) 



In view of (67), consider now a special case: s = 1/2. Then, the weight func- 
tion is known as lognormal distribution and polynomials Pn{u'. 1/2) (up to a 



Notice that a = I J dxq 
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constant ) arc known as Sticltjcs-Wigcrt (S-W) polynomials. Their physical rel- 
evance will be discussed below in Subsection 4.6. In the meantime, we introduce 
the Fourier transform of /(a;) in the usual way as 



oo 

/ 



dxexp{2Triex)f{x) = f{6) (68) 
Then, the Parseval relation implies: 

oo CO 

j ^n{x)^m{x)dx = J ^n{0)^m{0)dx = \\<^n{x)\f S,nn, (69a) 



causing 



^n{6) = (J) ' exp(- (tt/c) 0') C]:{q){-q-h'^'')K (69b) 



*:=0 



Comparison between these results and Eq.(32) produces 



CO 

/ 



Hr^i-q-h-^-''; q)Hrr,{-q-h-^^'';q) exp(-2 (tt/c) 0^) = (|-) ^ q-^{q, g)„(5. 



(70a) 



which can be alternatively rewritten as 
1 



J H^{-q-h-^^'')Hm{-q-h-^^'')M^Tre;q)de = q-''{q,q)nSmn (70b) 



with ^3{9,q) = q^ /2e*"^ . That is is one of the Jacobi's theta func- 



n= — OO 



tions. In order to use the obtained results, it is useful to compare them against 
those, known in literature already, e.g. see [54]. Our Eq.(70a) is in agreement 
with (5) of [54] if we make identifications: k = n and c = \/2«;, where k is the 
parameter introduced in this reference. With help of such an identification we 
can proceed with comparison. For this purpose, following [54] we introduce yet 
another generating function 



Sn{t;q) 



n 
k 



q" t" (71) 



fc=0 

so that the S-W polynomials can be written now as 

S^{t;q) = (-l)"g*(v/M^)-ip„(i;^) = {-iTq— {^/i^)-'Sr^{-qh■,q), 

(72) 
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in accord with Rcf. [55], page 197, provided that < g < 1. Comparison between 
generating functions (32) and (72) allows us to write as well 



Sn{t; q-^) = i?„(ig-"; q), or, equivalently, q-^) = 5„(Q-"i; q). (73) 

Using this result we can rewrite the recursion relation (33) for Hn{t; q) in terms 
of the recursion relation for Sn{t;q) if needed and then to repeat all the argu- 
ments with creation and annihilation operators, etc. For the sake of space, we 
leave this option as an exercise for our readers. Instead, to finish our discus- 
sion we would like to show how the obtained polynomials reduce to the usual 
Hermite polynomials in the limit — > 1~. For this purpose we would like to 
demonstrate that the recursion relation (33) is actually the recursion relation 
for the continuous g— Hermite polynomials [56,57]. This means that we have to 
demonstrate that under some conditions (to be specified) the recursion (33) is 
equivalent to 

2xHn{x I q) = Hn+i{x I ?) + (1 - q^Hn-iix I q). (74) 

known for q-Hermite polynomials. To demonstrate the equivalence, we assume 
that X = cos 6 in (74) and then, let ^ = e'^. Furthermore, we assume that 

Hn{x\q) = z^Hn{z-^;q) (75) 

allowing us to obtain 

{z + z-^)z^Hn = z^+^Hn+1 + (1 - (z")z"-^7f„_i. (76) 

Finally, we set = ^/t which brings us back to Eq.(33). This time, however, 
we can use results known in literature for g— Hermite polynomials [41-43] in 
order to obtain at once 



(77) 



where Hn{x) are the standard Hermitian polynomials. In view of Eq.s(72),(73), 
not surprisingly, the S-W polynomials are also reducible to Hn{x). Details can 
be found in the same references. 




4.4 ASEP, q-deformed hcirmonic oscillator and spin chains 

In this subsection we would like to connect the results obtained thus far with 
that for XXX and XXZ spin chains. Although a connection with XXX spin 
chain was established already at the beginning of this section, it is extremely 
helpful to arrive at the same conclusions using alternative (physically inspired) 
arguments and methods. To understand the logic of our arguments we encourage 
our readers to read Appendix at this point. In it we provide a self contained 
summary of results related to the asymmetric simple exclusion process (ASEP), 
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especially emphasizing its connection with static and dynamic properties of 
XXX and XXZ spin chains. 

ASEP was discussed in high energy physics literature, e.g. see Ref.[58], in 
connection with the random matrix ensembles. To avoid repeats, we would like 
to use the results of Appendix in order to consider a steady-state regime only. 
To be in accord with literature on ASEP, we complicate matters by imposing 
some nontrivial boundary conditions. 

In the steady -state regime Eq.(A.12) of Appendix acquires the form : SC = 
A. Explicitly, 

SC =plED- prDE. (78) 

In the steady-state regime, the operator S becomes an arbitrary c-number [59]. 
In view of this, following Sasamoto [60] we rewrite (78) as 

PrDE-plED = C{D + E). (79) 

Such operator equation should be supplemented by the boundary conditions 
which are chosen to be as 

a<W\E = C,<W\ and/3£»|y>=Cl^>. (80) 

The normalized steady-state probability for some configuration C can be written 
now as 

^^^^ <w\c^\v> ^^^^ 

with the operator Xi being either D or E depending on wether the i— th site is 
occupied or empty. To calculate P{C) we need to determine ^ while assuming 
parameters a and /3 to be assigned. In Appendix we demonstrate that it is 
possible to equate ( to one so that, in agreement with Ref.[61], we obtain the 
following representation for D and E operators: 

D = - + , ^ a, E 



1-q y/l-q l-q v^l - q 

converting equation (78) into (52d). In view of this mapping into g— deformed 
oscillator algebra, we can expand both vectors | V > and < W \ into Fourier 
series, e.g. | V >= J2rn^m{V) \ m > where, using (57), we put | m >= ^f^. 
By combining equations (57) and (80) and results of Appendix we obtain the 
following recurrence equation for f2„ : 

1 



^n{V){^ - 1) = nn+i{V)^l - (82) 

Following [47], we assume that < \ V >= 1. Then, the above recurrence 
produces 

V(9'9)" 
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with parameter v = ~ 1. Analogously, we obtain: 

nr,{w) = -==, 

with w = — 1. Obtained results exhibit apparently singular behavior for 
g — > 1~. These singularities are only apparent since they cancel out when one 
computes quantities of physical interest discussed in both [62] and [63]. As 
results of Appendix indicate, such a crossover is also nontrivial physically since it 
involves careful treatment of the transition from XXZ to XXX antiferromagnetic 
spin chains. Hence, the results obtained thus far enable us to connect the 
partition function (27) (or (32)) with either XXX or XXZ spin chains but are 
not yet sufficient for making an unambiguous choice between these two models. 
This task is accomplished below, in the rest of this section. 

4.5 Crossover between the XXZ and XXX spin chains: 
connections with the KPZ and EW equations and the 
lattice Liouville model 

Following Derrida and Malick [63] , we notice that ASEP is the lattice version of 
the famous Kradar-Parisi- Zhang (KPZ) equation [64]. The transition q ^ 1~ 
corresponds to a transition (in the sense of renormalization group analysis) from 
the regime of ballistic deposition whose growth is described by the KPZ equa- 
tion to another regime described by the Edwards- Wilkinson (EW) equation. 
In the context of ASEP (that is microscopically) such a transition is discussed 
in detail in [65]. Alternative treatment is given in [63]. The task of obtaining 
the KPZ or EW equations from those describing the ASEP is nontrivial and 
was accomplished only very recently [66,67]. It is essential for us that in doing 
so the rules of constructing the restricted solid-on -solid (RSOS) models were 
used. From the work by Huse [68] it is known that such models can be found in 
four thermodynamic regimes. The crossover from the regime III to regime IV 
is described by the critical exponents of the Friedan, Qui and Shenker unitary 
CFT series [69] . The crossover from the regime III to regime IV happens to be 
relevant to the crossover from KPZ to EW regime as we would like to explain 
now. 

As results of our Appendix indicate, the truly asymmetric simple exclusion 
process is associated with the XXZ model at the microscopic level and with 
the KPZ equation/model at the macroscopic level. Accordingly, the symmetric 
exclusion process is associated with the XXX model at the microscopic level 
and with the EW equation/model at the macroscopic level. At the level of the 
Bethe ansatz for open XXZ chain with boundaries full details of the crossover 
from the KPZ to EW regime were exhaustively worked out only very recently 
[70]. For the purposes of this work it is important to notice that for certain 
values of parameters the Hamiltonian of open XXZ spin chain model with 

l^That is Eq.(1.3) of [70]. 
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boundaries can be brought to the foUowing canonical form 
1 1 1 

(83) 

In the case of ASEP we put q — \J'PrIpl so that for physical reasons parameter q 
is not complex. However, mathematically, we can allow q to become complex. In 

particular, following Pasquer and Saleur [71], we can let q = e*'' with 7 = . 

/i + 1 

For such values of q use of finite scaling analysis applied to the spectrum of the 
above defined Hamiltonian produces the central charge 

c=l-^^,, = 2,3,.... (84) 

of the unitary CFT series. Furthermore, if is generator of the Temperley-Lieb 
algebr£0, then Hxxz can be rewritten as [72] 



Hxxz = -Y.^e^-\{q + q^')]. (85) 
i=i 

This fact allows us to make immediate connections with quantum groups and 
theory of knots and links. Below, in Section 5 we shall discuss briefly different 
arguments leading to analogous conclusions. 

The results just described allow us to connect CFT with the exactly inte- 
grable lattice models. In view of this, one can pose the following question: 
Given the connection we just described, can we write down explicitly the cor- 
responding path integral string-theoretic models reproducing results of exactly 
integrable lattice models at and away from criticality? Before providing the 
answer in the following subsection, we would like to conclude this subsection 
with a partial answer. In particular, we would like to mention the work by 
Faddeev and Tirkkonen [73] connecting the lattice Liouville model with the 
spin 1/2 XXZ chain. Based on this result, it should be clear that in the region 
c < 1 it is indeed possible by using the combinatorial analysis described above to 
make a link between the continuum and discrete Liouville theorieJ^. It can be 
made in such a way that, at least at crtiticality, the results of exactly integrable 
2 dimensional models are in agreement with those which are obtainable field- 
theoretically. The domain c > 1 is physically meaningless because the models 
{other than string-theoretic) we discussed in this section loose their physical 
meaning in this region. This conclusion will be further reinforced in the next 
subsection. 



4.6 ASEP, vicious random walkers and string models 

We have discussed at length the role of vicious random walkers in our treat- 
ment of the Kontsevich-Witten (K-W) model [74]. Forrester [75] noticed that 

^^That is = e^, eiei^iCi = q~^ei and eiej = e^ei for \i — j\ > 2. 
^^The matrix c = 1 theories will be discussed separately below. 
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the random turns vicious walkers model is just a special case of ASEP. Further 
details on connections between the ASEP, vicious walkers, KPZ and random 
matrix theory can be found in the paper by Sasamoto [76]. In the paper by 
Mukhi [77] it is emphasized that while the K-W model is the matrix model 
representing c < 1 bosonic string, the Penner matrix model with imaginary 
coupling constant is representing c = 1 bosonic Euclidean string on the cylinder 
of (self-dual) radius R ~ 10. Furthermore, Ghoshal and Vafa [79] have demon- 
strated that c = 1,R = 1 string is dual to the topological string on a conifold 
singularity. Wc shall briefly discuss this connection below. Before doing so, it 
is instructive to discuss the crossover from c = 1 to c < 1 string models in terms 
of vicious walkers. To do so we shall use some results from our work on K-W 
model and from the paper by Forrester [75] . 

Thus, we would like to consider a planar lattice where at the beginning 
we place only one directed path P: from (a, 1) to (6, Tvji^l- The information 
about this path can be encoded into multiset Hory{P) of y-coordinates of the 
horizontal steps of P. Let now 

wiP)= H X,. (86) 

i=Hory{P) 

Using these definitions, the extension of these results to an assembly of directed 
random vicious walkers is given as a product: W{P) = w(Pi) • ■ ■ w{Pk)- Finally, 
the generating function for an assembly of such walkers is given by 

p 

where W{P) is made of monomials of the type " provided that 

mi -|- • • • -f- TOjv — b ^ a. The following theorem discussed in our work, Ref[60], 
is of central importance for calculation of such defined generating function. 

Given integers < ai < ■ ■ ■ < Uk and < 5i < • • • < 6^ , let M^j be the 
k X k matrix M^j = hb^-ai{xi, ...,xn) then, 

detM = EM^(P), (88) 

p 

where the sum is taken over all sequences (Pi, ■■■,Pk) = P of nonintcrsecting 
lattice paths Pi : {ai, 1) — > {bi, N), i — 1 — k. 

Let now ai = i and bj = Xi -\- j so that 1 < i,j < k with A being a 
partition of TV with k parts then, detM = sx{xi, ...,xn), where sa(x) is the 
Schur polynomial. In our work [74] we demonstrated that in the limit iV — s- oo 
such defined Schur polynomial coincides with the partition (generating) function 
for the Kontsevich model. To obtain results of Forrester requires us to apply 
some additional efforts. These are worth discussing. 

^*This was initially demonstrated by Distlcr and Vafa [78]. 
^^Very much in the same way as wc have discussed in Section 2. 
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Unlike the K-W case, this time, we need to diseuss the continuous random 
walks in the plane. Let x-coordinate represent "space" while y-coordinate- 
"time". If initially {t = 0) we had fc- walkers in positions —L < x\ < X2 < 
■ ■ ■ < Xk < L, the same order should persist V t > 0. At each tick of the clock 
each walker is moving either to the right or to the left with equal probability p 
(that is we are in the regime appropriate for the XXX spin chain in the ASEP 
terminology). As before, let Xg — (xi,o, ...,Xkfi) be the initial configuration of 
A;— walkers and x/ = {xij, ...,Xkj) be the final configuration at time t. To 
calculate the total number of walks starting at time f = at xq and ending at 
time t at xy we need to know the probability distribution Wk (xq ; i) that 

the walkers proceed without bumping into each other. Should these random 
walks be totally uncorrelated, we would obtain for the probability distribution 
the standard Gaussian result 



exp{ -(x/-xo )V2£)t} 
(27rDi)* 



Wk (xo ^ X/, t) ^ fc/2 • (89) 



In the present case the walks arc restricted (correlated) so that the probability 
should be modified. This modification can be found in the work by Fisher and 
Huse [80]. These authors obtain 



^exp{-(x2 +xg)/2I?0 
(27r£)t) 

with 



Wki^o ^ xy;t) = C/fc(xo,xy;t) ;J^,^°/2 ' (9^) 



[/fe(xo,x/;i) = J2 e{g)eM^^^^^p^]- (91) 

geSk 

In this expression e{g) = ±1, and the index g runs over all members of the sym- 
metric group Sk- Mathematically, following Gaudin [81], this problem can be 
looked upon as a problem of a random walk inside the A:— dimensional kaleido- 
scope (Weyl cone) usually complicated by imposition of some boundary condi- 
tions on walks at the walls of the cone. Connection of such random walk problem 
with random matrices was discussed by Grabiner [82] whose results were very 
recently improved and generalized by Krattenthaller [83]. In the work by de 
Haro some applications of Grabiner's results to high energy physics were con- 
sidered [84] . Here we would like to approach the same class of problems based 
on the results obtained in this paper. In particular, some calculations made in 
[80] indicate that for L — > 00 with accuracy up to O (L'^/Dt) it is possible to 
rewrite ?7fe (xq , x/ ; t) as follows: 

?7fc(xo,x/;t) ~ consiA(x/)A(xo)/ {Dt)'"' + 0{L'^/Dt) (92) 

with uk = (1/2) k{k - 1) and const = 1/1!2! ■ ■ ■ {k - 1)1 and A(x) being the 
Vandermonde determinant, i.e. 

A(x) = Hix, - xj). (93) 

i<j 



24 



Next, from standard texts in probability theory it is known that non-normalized 
expression for VF°(xo ^ x/; in the Umit of long times provides the number of 
random walks of n steps (since n ^ t) from point xq to point x/. Hence, the 
same must be true for Wki^o ^/S ^-nd, therefore, Wk{xo — > x/; t) « det M. 
Consider such walks for which x{t = 0) = Xq = x.{t = tf) = x/. Then, using 
(90) and (92) we obtain the probability distribution for the Gaussian unitary 
ensemble [85], i.e. 

Wk{-^o = X/; t) = const A2(x) exp(-x2). (94) 

Some additional manipulations (described in our work [74]) using this ensem- 
ble lead directly to the K-W matrix model. Forrester [75] considered a related 
quantity- the probability that all k vicious walkers will survive at time tf. To 
obtain this probability requires integration of Wk (xq — + ^f]t) over the simplex 
A defined hy — L < xi < X2 < ■ ■ ■ < Xk < lF°l. Without loss of generality, it is 
permissible to use Wk (xq = Xf;t) instead of Wk (xq ^f',t) in calculating such 
a probability. Then, the obtained result coincides (up to a constant) with the 
partition function for topological gravity Z{g), Eq.(3.1) of [86.87p^. Further- 
more, such defined partition function can be employed to reproduce back the 
Hermite polynomial Hk{x) defined by (77) which has an interpretation as the 
wavefunction(amplitudc) of the FZZT Z?— brane [86,87]. Specifically, we have 

< det(x-A/) (|)^i7fc(a;yi) 

= ^ j dM det(x - M)e- s*"^*^' . (95) 

This expression is a special case of Heine's formula representing monic orthog- 
onal polynomials through random matrices. In the above formula k is related 
to the size of Hcrmitian matrix M and g is the coupling constant. 

Following Forrester [75], the result, Eq.(90), can be treated more accurately 
(albeit a bit speculatively) if, in addition to the parameter D we introduce 
another parameter a - the spacing between random walkers at time t — Q. 
Furthermore, if the time direction is treated as space direction (as it is commonly 
done for Id quantum systems in connection with 2d classical systems), then yet 
another parameter T{k,t) should be introduced which effectively renormalizes 
D. This eventually causes us to replace Z{g) by the following integral (up to a 
constant) 

fe °° 

i(<?) J dx, exp(-^ In^ X,) l[{x, - x^f = j dMe-^'^^"'^'^\ (96) 



^"Such type of integration is described in detail in our papers, Parts I and II, from which 
it follows that in the limit L — » oo such a simplex integration can be replaced by the usual 

fc oo k 

integration, i.e Yi ' " <^ H / 11 ' " accord with Forrester. 

i —ooi 

■^^Since the hermitian matrix model given by Eq.(3.1) is just a partition function for the 
Gaussian unitary ensemble [85]. 
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Tierz demonstrated [88] that Z{g) (up to a constant) is partition function of 
the Chern-Simons (C-S) field theory with gauge group U{k) hving on 3-sphere 
S^. Okuyama [87] used (95) in order to get analogous result for a D-brane 
amplitude in the C-S model. Using Heine's formula, he obtained the Stieltjes- 
Wiegert (S-W) polynomial, our Eq.(71), which can be expressed via Rogers- 
Szego polynomial according to (73) and, hence, via the g— Hermite polynomial 
in view of the relation (75). Since in the limit g — > 1~ the 5— Hermite polynomial 
is reducible to the usual Hermite polynomial according to (77) , there should be 
analogous procedure in going from the partition function Z(g) to Z{g). Such a 
procedure can be developed, in principle, by reversing arguments of Forrester. 
However, these arguments are much less rigorous and physically transparent 
than those used in previous subsection where we discussed the crossover from 
XXZ to XXX model. In view of the results presented below, in the following 
section, we leave the problem of crossover between the matrix ensembles outside 
the scope of this paper. To avoid duplications, we refer our readers to the 
paper by Okuyama [87] where details are provided relating our results to the 
topological A and B -branes. 

5 Discussion. Gaudin model as linkage between 
the WZNW model and K-Z equations 

We would like to remind to our readers that all results obtained thus far 
can be traced back to our Eq.(32) defining the Rogers-Szego polynomial which 
physically was interpreted as partition function for the ferromagnetic P-F spin 
chairF^. In previous section numerous attempts were made to connect this par- 
tition function to various known string models, even though already in Section 
4.1 we came to the conclusion that in the limit of infinitely long chains the 
antiferromagnetic version of P-F spin chain can be replaced by the spin 1/2 
antiferromagnetic XXX chain. If this is so, then from literature it is known 
that the behavior of such spin chain is described by the SUi{2) WZNW model 
[44]. Hence, at the physical level of rigor the problem of connecting Veneziano 
amplitudes to physical models can be considered as completely solved. Here 
we argue that at the mathematical level of rigor this is not quite so yet. This 
conclusion concerns not only problems discussed in this paper but, in general, 
the connection between WZNW models, spin chains and K-Z equations. 

It is true that the K-Z equations and WZNW model are inseparable from 
each other [44] but the extent to which spin chains can be directly linked to 
both the WZNW models and K-Z equations remains to be explained. Following 
Varchenko [89], we notice that the link between the K-Z equations and WZNW 
models can be made only with help of the Gaudin model, while the connection 
with spin chains can be made only by using the quantum version of the K- 

^^The antiferromagnetic version of P-F spin chain is easily obtainable from this ferromag- 
netic version as discussed in Section 3. 
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Z equation. Such quantized version of the K-Z equation is not immediately 
connected with the standard WZNW model as discussed in many places [89, 
90] . For the sake of space we refer our readers to these references and to our 
latest work connecting Vcneziano amplitudes with Gaudin models [27] in which 
we also discuss in detail physical implications of such a connection. 



Appendix. Basics of ASEP 

A.l. Equations of motion and spin chains 

The one dimensional asymmetric simple exclusion process (ASEP) had been 
studied for some time [91]. The purpose of this Appendix is to summarize the 
key features of this process which are of immediate relevance to the content of 
this paper. To this purpose, following Schiitz [92], we shall briefly describe the 
ASEP with sequential updating. Let Bn ■= {xi, xn} be the set of sites of one 
dimensional lattice arranged at time t in such a way that xi < X2 < ■ ■ ■ < Xn- 
It is expected that each time update will not destroy this order. 

Consider first the simplest case oiN = 1. Let pn {p^) be the probability of 
a particle located at the site x to move to the right (left) then, after transition to 
continuous time, the master equation for the probability P(x; t) can be written 
as follows 

^P(a;; t) = prP{x - l;t) + plP{x +l;t)- P{x; t). (A.l) 
Assuming that P{x; t) = exp{—et)P{x) so that that 

27r 

P{x;t) = J dpexp{—et)f{p)exp{ipx), p G [0,27r), we obtain the dispersion 



relation for the energy e{p) : 

e{p)=pR{l-e-'n+PLil-e'n- (A.2) 

The initial condition P{x;0) = Sx,y determines the amplitude f{p) = e~'P^/27r 
and finally yields 

P{x; t\y;0) = ^ [ rfpe-^^^ )*e-'f^e'f^ = e-(9+«"')^*g^-^4_j,(2£)i), (A.3) 
2tt J 



where q = \/pr/pl, D = ^prPl and /„(2Dt) is the modified Bcssel function. 
These results can be easily extended to the case N = 2. Indeed, for this case 
we obtain the following equation of motion 

£P{Xi,X2) = -pr{P{xi-1,X2)+P{xi,X2-1)-2P{xi,X2)) 

-Pl{P{xi + 1, X2) + P{xi,X2 + 1) - 2P(xi, a;2)) (A.4) 
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which should be supplemented by the boundary condition 



P{x, a; + 1) = prP{x, x) + plP{x + 1, a; + 1) Vx. (A.5) 

Imposition of this boundary condition allows us to look for a solution of (A. 4) 
in the (Bethe ansatz) form 

Pixi,X2) = ^lae^Pi^^e'f^^^ + ylaie^P^^^e'^ (A.6) 

yielding e{pi,p2) ~ ^{pi) + e(P2)- Use of the boundary condition (A.5) fixes the 
ratio (the S-matrix) 6*12 = ^12/^21 as follows: 

S{pi,P2)^ ■ -— Z^; — . A.7) 

PR + ple P^ P^ ^ 

To connect this result with the quantum spin chains, we consider the case of 
symmetric hopping first. In this case we have pn = pl — 1/2 so that (A.7) is 
reduced to 

1 + e*Pi+*P2 _ 2e^P^ 
5'xxx(pi,P2) = -Y:j7^7^7T^^n7^;^ (A.8) 

from this result we can recognize the S matrix for the XXX spin 1/2 Heisenberg 
ferromagnet [81]. If p_R 7^ pL, to bring (A.7) in correspondence with the spin 
chain S— matrix requires additional efforts. Following Gwa and Spohn [93] we 
replace the complex numbers e*^^ and e'P^ in (A.7) respectively by zi and Z2. 
In such a form (A.7) exactly coincides with the S-matrix obtained by Gwa 

and SpohrF^. After this we can rescale z^ {i = 1,2) as follows: Zi ~ \J^^i- 

Substitution of such an ansatz into (A.7) leads to the result 

Sxxz{zi,z2) = -7—^-: ^-r-^, (A.9) 

1 + Z1Z2 — 2.i\Z2 

provided that 2A = ^1 x^/vlPr- For pR = pl — 1/2 we obtain A = 1 as 
required for the XXX chairo Ifi however, pn 7^ p^, then, the obtained S- 
matrix coincides with that known for the XXZ model [81] if we again relabel z^ 
by e^P^ which is always permissible since the parameter p is determined by the 
Bethe equations (to be discussed below) anyway. 

In the case of XXZ spin chain it is customary to think about the massless 
— 1 < A < 1 and massive |A| > 1 regime. The massless regime describes various 
CFT discussed in the text while the massive regime describes massive excitations 
away from criticality. As Gaudin had demonstrated [81], for XXZ chain it is 
sufficient to consider only A > domain which makes XXZ model perfect for 
uses in ASEP. The cases A = and A — > 00 also physically interesting: the 
first corresponds to the XY model and the second to the Ising model. 



23E.g. see their Eq.(3.5). 

^^It should be noticed though that such a parametrization is not unique. For instance, 
following [70] it is possible to choose a slightly different parametrization, e.g. A = — ^((j+q"^), 
where q = ^Pr/pl- 



28 



Once the S-matrix is found, the TV— particle solution can be easily con- 
structed [92]. For instance, for = 3 we write 

\l/(a;i, X2, xs) = exp(ipia;i + 1^22:2 + ipsxs) + S21 e^p{ip2Xi + ipiX2 + ips^s) 

+'S'32S'3i exp(ip2a;i + 1^32:2 + ipixs) + S21S31S32 exp{ip3Xi + ip2X2 + ipsxs) 
+S31S32 ex.p{ip3Xi + ipiX2 + ip2X3) + S32 exp{ipiXi + ip3X2 + ip2X3), 

etc. This result is used instead of f{p) in (A. 3) so that the full solution is given 
by 

P{xi, xn; t I 2/1, 0) = TT / dpze-^'P'^'e-'P'^'^-lxi, xn)- 

(A.ll) 

The above picture can be refined as follows. First, the particle sitting at Xi will 
move to the right(left) only if the nearby site is not occupied. Hence, the prob- 
abilities PR and PL can have values ranging from to 1. For instance, for the 
totally asymmetric exclusion process (TASEP) particle can move to the right 
with probability 1 if the neighboring site to its right is empty. Otherwise the 
move is rejected. Since under such circumstances particle can never move to 
the left, there must be a particle source located next to the leftmost particle 
position and the particle sink located immediately after the rightmost position. 
After imposition of emission and absorption rates for these sources and sinks, 
we end up with the Bethc ansatz complicated by the imposed boundary condi- 
tions. Although in the case of solid state physics these conditions are normally 
assumed to be periodic, in the present case, they should be chosen among the 
solutions of the Sklyanin boundary equations [59, 72]. At more intuitive level of 
presentation compatible with results just discussed, the Bethe ansatz for XXZ 
chain accounting for the boundary effects is given in the pedagogically written 
paper by Alcaraz et al [94]. 



A. 2. Dynamics of ASEP and operator algebra 

To make these results useful for the main text, several additional steps are 
needed. For this purpose we shall follow works by Sasamoto and Wadati [95] 
and Stinchcombe and Schiitz [59] . In doing so we rederive many of their results 
differently. 

We begin with observation that the state of one dimensional lattice contain- 
ing A'^ sites can be described in terms of a string of operators D and E, where 
D stands for the occupied and E for empty fc— th position along Id lattice. The 
non normalized probability (of the type given in (A.ll)) can then be presented 
as a sum of terms like this < EEDEDDD ■ ■ ■ E > to be discussed in more 
detail below. 

Let C = D + E he the time- independent operator. Then, for the operator 
D to be time-dependent the following commutation relations should hold 

SC + DC = A, (A.12a) 
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CS-CD = A, (A.12b) 

DD + Di) =[D,S]. (A.12c) 

If A = plCD -prDC +{pr - pl)D'^ (or A = plED- prDE ), it is possible 
to determine S using equations (A. 12) so that we obtain, 

b=\[K,C-\ (A.13a) 

5-i{A,C-i}, (A.iab) 

provided that 

KC-^D ^ DC-^A (A.13c) 

with { , } being an anticommutator. 

As before, we want to consider the case pr — pj^ =i first. This condition 
leads to A = i [C, D] . It is convenient at this stage to introduce an operator D„ 
^fjn-i ]j(j-ri ^j^^ Fourier transform Dp —J2n exp{ipn). Using (A. 13a) 
with A just defined leads to the following equation of motion for £)„ : 

I)„ = ^[Z?„+i+D„_i-2D„] (A.14) 

to be compared with (A.l). Such a compariso n p roduces at once 'Dp{t) = 
exp(— e(p))I?p(0) so that e{p) — 1 — cos v as beforj^. Consider now Eq. (A. 13c). 
Under conditions p_r = p_l =^ it can be written as 
CDC-^D +DC-^DC =2L>2 or, equivalently, a@ 

A. + DnDn-i = 2D,,D^. (A.15) 

Following Gaudin [81], we consider a formal expansion DnDm =aexp(ipi7i + 
ip2'm) + (3 eyip{ip2n + ipim) and use it in the previous equation in order to 
obtain: 

aexp(ipi(n + 1) + ip2n) + (3 eicp{ip2{n + 1) + ipin) 
+Q;exp(ipin + ip2{n — 1)) + /3 exp{ip2n + ipi{n — 1)) 
= 2aeiq){ipin + ip2n) + 2f3 exp{ip2n + ipin) . (A. 16) 

From here we also obtain: 



{a e-xjp{ipin+ip2n){eiip{ipi)+exp{—ip2)~2)+P exp{ipin+ip2n){exp{ip2)+exp(—ipi)—2) — 
and, therefore, 

a, ^ _ 1 + cxp(z(pi +P2)) ~2exp{ipi) /A1-7^ 

S{pi,p2) = -3 = 7T7 ■ TT — ^ — — exp I p2 -Pi (A.17) 

[3 1 + exp(i(pi + P2)) - 2 exp{ip2) 

25E.g. see (A.2) with = = 1/2 

^^Since using definition of Z?„ we have: -D„-|_i = CD„C~^ and D„_i = D^C. 
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to be compared with (A. 8). An extra factor exp(i(p2 — Pi)) can be actually 
dropped from the S*— matrix in view of the following chain of arguments. 
Introduce the correlation function as follows 

P{xi,...,XN;t I yi,...,yN;0) = Z]:,^Tr[Di{t)---DN{t)C^] 

N 2^ 

= 11^ I dpie-<P^')'e-'P^y^'i>{pi,...,PN), (A.18) 
1=1 

where ^{pi,...,pn) = Z]:^^Tr[Dp,{0) ■ ■ ■ Vp„{t)C^] and Zn = tr[C^]. In 
arriving at this result the definition of 'Dp{t), was used along with the fact that 
CVpC^^ ~ e^^'f Vp. Also, the invariance of the trace under cyclic permuta- 
tions and the translational invariance of the correlation function implying that 
Pn) 7^ only if pi = was taken into account. These conditions 
are sufficient for obtaining the Bethe ansatz equations 

N 

expfeTV) ^ J] Vi ^ j, (A.19) 

where S{pi,pj) is the same S— matrix as in (A. 17), except of the missing factor 
exp(j(pi — Pj)) which is dropped in view of translational invariance^. 

Extension of these results to the case pa ^ pl is nontrivial. Because of 
this, we would like to provide some details not shown in the cited references. In 
particular, contrary to claims made in [95], we would like to demonstrate that 
the system of equations (A. 12) obtained in [59] is equivalent to the system of 
equations 

[C, S] = 0, (A.20a) 

CD + CT - SD = -plCD +prDC, (A.20b) 

DC + DS -TC ^ PlCD -prDC, (A.20c) 

DD + DD =[T,D] (A.20d) 

obtained in [95] with the purpose of describing asymmetric processes. 

To make a comparison between (A. 12) and (A. 20), we notice that (A. 20) has 
operators S and T which cannot be trivially identified with those present in 
(A. 12). Hence, the task lies in making such an identification. For this purpose, 
if we assume that S in (A. 20) is the same as in (A. 12) then, in view of (A. 20a), 
by subtracting (A. 12b) from (A. 12a) we obtain: 

DC + CD = 0. (A.21) 

This leads to either D = C'^DC or t) = -CDC''^. Therefore, taking into ac- 
count that, by construction, C is time-independent, we obtain: D — —CDC~^ + 
O, where Q is some diagonal time-independent matrix operator. 



^^Surely, in the case when the effects of boundaries should be accounted, this factor should 
be treated depending on the kind of boundary conditions imposed. 
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Next, using these results we multiply (A. 20b) from the right by and 
(A. 20c) by from the left, and add them together in order to arrive at 
equation (19) of [95], i.e. 

2D = prC-^DC + plCDC-^ - {pR + Pl)D. (A.22) 

Also, by multiplying this result from the right by D we obtain equation (20) of 
[95], that is 

= prDC-^DC + plCDC-^D - {pr+pl)D^, (A.23) 

provided that [T, D] = 0. That this is indeed the case can be seen from the 
same reference where the following result for T is obtained: 

2T ={2+pr-pl)D+prC-^DC -plCDC~\ (A.24) 

Using it, we obtain: [T, £>] = 0, in view of the fact that [C^^DC, D] = and 
[CDC^^^ , D] — since D — ~CDC^^ + 6 as we have already demonstrated. 
Furthermore, (A.24) can be straightforwardly obtained by subtracting (A. 20c) 
(multiplied by from the right) from (A. 20a) (multiplied by C^^ from the 
left). Thus, contrary to the claims made in [95], equations (A. 12) and (A. 20) 
are, in fact, equivalent. Nevertheless, as claimed in [95], the system of equations 
(A. 20) is easier to connect with the Bethe ansatz formalism. 

Indeed, using already known fact that _D„ —C^^^DC^^^ Eq.(A.22) can be 
brought into the form: 

Dn = ^[pRDn+l + PLDn-l - iPR+PL)D„]. (A.25) 

This result is formally in agreement with previously obtained (A. 14) for the 
fully symmetric case. The authors of [95] have chosen such a normalization 
for probabilities pr and pl that for symmetric case pr — pl = 1 (instead of 
PR = PL = 1/2). To restore the normally accepted condition pR = pi^ = 1/2 
requires only to rescale time appropriately. This observation is consistent with 
the fact that the analog of equation (A. 15) (which plays the central role in Bethe 
ansatz-type calculations) obtained with help of (A.23) is given by 

PL Ai+l Ai + PRDnDn-l = ipR + PL)DnDn (A. 26) 

which holds true wether we choose Pr = Pl = ^ or pr = p^ = 1/ These 
results allow us to reobtain the 5— matrix for the XXZ model in a way already 
described. 

A. 3. Steady- state and q- algebra for the deformed harmonic os- 
cillator 

^*It should be noted though that the authors of [95] have erroneously obtained (e.g. see 
their equation (23)) pjj-D^ +Pl^j5+i = {pr + VL)DnDnj^i instead of our (A. 26). 
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Using (79) we write 



PrDE-plED^C{D + E). (A.27) 

Let now D — Ai + Bia and E = A2 + -620^ where Ai and Bi, i = 1,2, are 
some c— numbers. Substituting these expressions back to (A.27) we obtain the 
following set of equations 

C{Ai+A2)-eAiA2^C, (A.28a) 
where C is some constant to be determined below, and 

(Bi = eBiA2, (A.28b) 

CB2=eB2Ai. (A.28c) 

From here we obtain: Ai = A2 = A = C/e, with Bi, B2 being yet arbitrary 
c- numbers and s — pr ~ pl- We can determine these numbers by comparing 
our results with those in [61]. This allows us to select Bi = B2 = ^i-q ' 

— = — ~ C. q ~ — so that we obtain 

e ^-'i PR 

D = + Q, (A.29a) 



l-q 
1 1 



E = + Q+ (A.29b) 



and, finally, 

aa+ - ga+a = 1 (A. 29c) 

in accord with (52d). 
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